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We introduce the quantum stochastic walk (QSW), which determines the evolution of generalized 
quantum mechanical walk on a graph that obeys a quantum stochastic equation of motion. Using an 
axiomatic approach, we specify the rules for all possible quantum, classical and quantum-stochastic 
transitions from a vertex as defined by its connectivity. We show how the family of possible QSWs 
encompasses both the classical random walk (CRW) and the quantum walk (QW) as special cases, 
but also includes more general probability distributions. As an example, we study the QSW on a 
line, the QW to CRW transition and transitions to genearlized QSWs that go beyond the CRW and 
QW. QSWs provide a new framework to the study of quantum algorithms as well as of quantum 
walks with environmental effects. 



Many classical algorithms, such as most Markov-chain 
Monte Carlo algorithms, are based on classical random 
walks (CRW), a probabilistic motion through the vertices 
of a graph. The quantum walk (QW) model is a unitary 
analogue of the CRW that is generally used to study and 
develop quantum algorithms [1, 2 , 3J. The quantum me- 
chanical nature of the QW yields different distributions 
for the position of the walker, as a QW allows for super- 
position and interference effects |3]. Algorithms based 
on QWs exhibit an exponential speedup over their clas- 
sical counterparts have been developed 13 |6j [7]. QWs 
have inspired the development of an intuitive approach to 
quantum algorithm design [i8j, some based on scattering 
theory [3] . They have recently been shown to be capable 
of performing universal quantum computation |10j . 

The transition from the QW into the classical regime 
has been studied by introducing decoherence to specific 
models of the discrete-time QW [III [HI [13l [M]. De- 
coherence has also been been studied as non-unitary ef- 
fects on continuous-time QW in the context of quantum 
transport, such as environmentally-assisted energy trans- 
fer in photosynthetic complexes [T5J Uni [T7J [T5J [T^ and 
state transfer in superconducting qubits [201 121] ■ For 
the purposes of experimental implementation, the ver- 
tices of the graph in a walk can be implemented using 
a qubit per vertex (an inefficient or unary mapping) or 
by employing a quantum state per vertex (the binary 
or efficient mapping). The choice of mapping impacts 
the simulation efficiency and their robustness under de- 
coherence [22l [23l [21] • The previous proposed approaches 
for exploring decoherence in quantum walks have added 
environmental-effects to a QW based on computational 
or physical models such as pure dephasing [17j but have 
not considered walks where the environmental effects are 
constructed axiomatically from the underlying graph. 

In this work, we define the quantum stochastic walk 
(QSW) using a set of axioms that incorporate unitary 
and non-unitary effects. A CRW is a type of classical 
stochastic processes. From the point of view of the theory 



of open quantum systems, the generalization of a classical 
stochastic process to the quantum regime is known to 
be a quantum stochastic process pi. '251, E71 Egl 
which is the most general type of evolution of a density 
matrix, not simply the Hamiltonian process proposed by 
the QW approach. The main goal of this paper is to 
introduce a set of axioms that allow for the construction 
of a quantum stochastic process constrained by a graph. 
We call all the walks that follow these axioms QSWs. 
We will show that the family of QSWs includes both the 
CRW and the QW as limiting cases. The QSW can yield 
new distributions that are not found either in the CRW 
or the QW. The connection between the three types of 
walks discussed in this manuscript is summarized in Fig. 
[l] For clarity, we focus on continuous-time walks, but 




FIG. 1: The quantum stochastic walk (QSW) is a quantum 
stochastic process defined by a set of axioms that describe 
the connectivity of an underlying graph. The family of pos- 
sible QSWs encompasses the classical random walk (CRW) 
and the quantum walk (QW) as limiting cases. QSWs can be 
used to generate walks that continuously interpolate between 
QW and CRW, and therefore can be employed to analyze 
algorithmic classical-to-quantum transitions. The family of 
QSWs also describes more general quantum stochastic pro- 
cesses with a behavior that cannot be recovered by the QW 
and CRW approaches. 

also sketch the corresponding procedure for the discrete- 
time walks. The QSW provides the fundamental tools 
to study the quantum-to-classical transition of walks, as 
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well as new methods of control by the application of non- 
unitary operations on a quantum system. 

The CRW describes probabilistic motion over a graph. 
The dynamics of the probability distribution function is 
given by the master equation for a classical stochastic 
process, 

|p„ = 5]M,>fc, (1) 

b 

where the vector element ph is the probability of being 
found at the vertex b of the graph. The matrix M is the 
generator for the evolution; its structure is constrained by 
axioms derived from the connectivity of the graph. For 
example, if vertices a and b are connected, = —7, if 
they are not, — 0, and = da'y where da is the 
degree of vertex a. 

In analogy to the CRW, the QW [T] has been defined 
so that the probability vector element Pa is replaced with 
(a|'0), which evolves according to the Schrodinger equa- 
tion, 

^^{a\,p)^-zJ2{a\H\b)m, (2) 

b 

where is a Hamiltonian to be defined based on ax- 
ioms coming from the graph. A choice of this definition 
is {a\H\b) ~ . This unitary evolution effectively ro- 
tates populations into coherences and back [35]. The 
QW fails to completely capture the stochastic nature of 
CRW. The random aspect of the QW comes solely from 
the probabilistic nature of measurements performed on 
the wavefunction. 

Since a classical stochastic process can be generalized 
to the quantum regime by means of a quantum stochastic 
process, a CRW should be generalized to a QSW derived 
from the graph |33j . For the generalization, we identify 
the probability vector with elements Pa with a density 
matrix with elements paa, and generalize the evolution 
to a quantum stochastic process, — M[p\, where M 
is a superoperator [25l [26l |28] . To make this evolution 
look similar to Eq. ([ij, we write the density matrix in 
terms of its indices, p = X^a a ^^'^ quan- 
tum stochastic master equation becomes, 

with the tensor A^g^ = {a\M[ |fe)(/3| ] |a). This cor- 
respondence was pointed out by Mohseni et al. in the 
context of energy transfer |16j . 

For a quantum stochastic process to be related to a 
walk, the superoperator M. must reflect the graph. The 
connectivity of the vertices will impose conditions on the 
transition rates of M. Since the quantum stochastic 
process is more general than both the classical stochas- 
tic process and the Schrodinger equation, the correspon- 



dence of the connectivity of the graph to the rules im- 
posed on M should include and go beyond the connec- 
tivity axioms for each of those. For a vertex m connected 
to vertices that include vertex n, we define a processes 
|m) ^ \n) which occurs at some rate that can evolve |m) 
to and from \n) . Transition rates for vertices that are 
not connected are defined to be zero. We employ these 
connectivity rules as the main principle for defining valid 
QSWs from a given graph. To further explore the con- 
nection from the QSW to the CRW and QW as well as 
more general behaviors, we discuss the different limiting 
cases. 

For the classical case, the allowed transitions come 
from incoherent population hopping of the form 
|m)(r7i| ^ |7i)(n|, and, for completeness, |m)(m| ^ 
\m){m\. These conditions constrain M. to operate only 
on the diagonal elements of p, like Eq. ([T]). In other 
words, the QSW transition tensor must have the prop- 

The QSW should also recover the QW, where evolution 
among vertices happens through coherences developed 
by a Hamiltonian. These transitions include terms of the 
form, |m)(r7i| ^ |m)(n|, that exchange populations to co- 
herences with the connected vertices. For completeness, 
we also consider transitions of the form |m)(ri| ^ |m)(n|. 
If m is also connected to another vertex I, additional 
transitions can happen. One transition exchanges popu- 
lations into coherences among the two vertices connected 
to m, |rn)(m| ^ 10 ("•li the other exchanges coherences 
between m and a connected vertex into coherence among 
two vertices connected to m, |m)(n| ^ lOi*^!- These con- 
ditions allow for the recovery of Eq. Q. By including 
the conjugates of these, we have now exhausted all the 
possibilities of the transitions that can happen following 
the connectivity of vertex m. These rules can be applied 
to vertices with any number of connections, and serve 
as the basis for the correspondence of the graph to the 
QSW. 

To complete the generalization, we need an equation 
of motion for a quantum stochastic process. We choose 
the Kossakowski-Lindblad master equation |26j [27l [28] . 
ip_^ C[p\ = Y.k-m^P]-\L^uLkP-\pLlU + LupLl 
which evolves the density matrix as a quantum stochas- 
tic process inspired by environmental effects under the 
Markov approximation. The Hamiltonian term describes 
the coherent evolution (Schrodinger equation) while the 
rest describe a stochastic evolution. We now set M ^ C 
in Eq. (|3|, where, 

+ Sab{f3\ (^iH-^LlLk^ \a) + {a\Lk\b){l3\Ll\a). 

The connectivity conditions between a vertex m and 
some connected vertices n and I and the corresponding 
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TABLE I: Axioms for the quantum stochastic walk for processes that connect vertex \m) to its neighbors. Sum over k, and 
conjugate elements are implied. Axioms (1) and (2) correspond to the classical random walk. Axioms (3), (4) and (5) contain 
the quantum walk using H and additional evolution due to {Lk} terms. Axiom (6) comes only from the quantum stochastic 
walk, having no equivalent in the classical random walk or quantum walk. 



non-zero transition rates according to Eq. Q can be sum- 
marized in Table |T] pH;. 

The Axioms from Table [T] capture the behavior 
sketched in Fig. ([ij. To recover the CRW, it suffices to 
consider Axioms (1) and (2). These Axioms are a clas- 
sical subset of the transition rates of the QSW. On the 
other hand, the QW is obtained by making all the rates 
of each element of {Lk} zero. In this case, only the subset 
of Axioms (3), (4) and (5) from Table |l] are relevant, cor- 
responding to the Hamiltonian, and the QW is recovered. 
If the rates of {Lk} are nonzero, these Axioms contain 
behavior beyond the QW. Finally, Axiom (6) has prop- 
erties that have no equivalent in either the CRW or the 
QW; it is a type of transition that appears exclusively in 
the QSW and leads to different distributions. The choice 
of H or {Lk} is not uniquely determined by the Axioms. 

For example, a CRW is equivalent to a QSW with the 
choice of each element of {Lk} to be associated to M in 
the following manner. First, since connections are de- 
fined between two vertices, it is easier to write the index 
k in terms of two indices k = (k, k'). Using this notation, 
we enforce the graph by choosing ^(k.k') = \k){k'\. 
The Hamiltonian rates of H are set to zero. This choice 
ensures all the transition rates to be zero except Axiom 
(2) from Table |l] thereby recovering the CRW. An inter- 
polation between the CRW and the QW conditions can 
be used to obtain the classical-to-quantum transition of 
any graph by changing the relative rates of the unitary 
and non-unitary processes. 

Other sets of {Lk} can yield behavior different from 
both the classical and quantum walks. To illustrate the 
difference between these choices, we study a simple ex- 
ample: the walk on an infinite line, where a vertex j is 
only connected to its nearest neighbors. The conditions 
for the dynamics of the walk on the line can be obtained 
from Table [l] by making j — m, j — 1 = I and j + 1 — n. 
To obtain behavior that is different from both the CRW 
and the QW, we specify a set {Lk} with only one member 
that is L = k' This choice makes Axiom 

(6) nonzero, a type of transition that cannot be inter- 



preted either as the CRW or the QW. We illustrate the 
resulting distribution, and compare it to CRW and the 
QW, in Fig.[2| 

The QSW provides a general framework to study the 
transition between the CRW and the QW. In Fig. |2^, 
we use the parameter lu = [0, 1] to interpolate between 
the QW and the CRW: C^[p] = -(1 - uj)i[H,p] + 

^ J2k i^^\^\^kP - \pL\Lk + LkpLl^ . By combining 
the unitary evolution and the evolution due to the 
{L{k,k')} terms, the transition from the QW to the CRW 
can be studied as one quantum stochastic walk. This 
procedure can be done in general to study the transition 
from any Q W to a CRW for any graph [35j . To highlight 
the difference between the QSW and the QW, we show 
the transition between them with L ~ k' ^'^k 
on Fig. [2]d. 

Although in this paper we focused on continuous-time 
walks, a parallel argument holds for discrete-time walks. 
A discrete-time CRW evolves each time step by a Markov 
chain with a stochastic matrix S following pa = pb- 
The quantum analogue is to use a quantum stochastic 
map S [25] by means of p'^^ = J^tp^bp PbfS, or equiv- 
alently as a superoperator, p' = B[p] = J^k^kpCj.- A 
similar set of Axioms to Table |l] can be computed by 
using B^^ = Efe(a|Cfe|&)(/9|C^|a) The connection 

between the discrete-time QW and the continuous-time 
Q W has been studied by Strauch [5IT| . 

In conclusion, we introduced an axiomatic approach to 
define the quantum stochastic walk, which behaves ac- 
cording to a quantum stochastic process as constrained 
by a particular graph. This walk recovers the classical 
random walk and the quantum walk as special cases, 
as well as the transition between them. The quantum 
stochastic walk allows for the construction of new types 
of walks that result in different probability distributions. 
As an example, we studied the walk on a line. The quan- 
tum stochastic walk provides a framework to study quan- 
tum walks under decoherence. Reexamination of pre- 
vious work that considered environmental effects from 
physical motivations might suggest that, if interpreted 



4 



0.25 



0.5 
5 15- <e< 




0.05- 



10 

Position 



FIG. 2: Probability distributions of a walker starting at po- 
sition of an infinite line defined by M are shown at time 
t = 5. The heavy line represents the quantum walk, the thin 
line represents the classical random walk and the dotted line 
represents a choice of a quantum stochastic walk due to a 
single operator L = ^, whose distribution can- 

not be obtained from the classical or quantum walks, a) The 
transition of the quantum walk to the classical random walk 
is shown. The diagonal elements of the density matrix (Pop- 
ulation) are plotted as a function of the vertex of the walker 
(Position) and a coupling parameter uo. When = the evo- 
lution is due only to the Hamiltonian, and when u = 1 the 
Hamiltonian dynamics are not present with the environmen- 
tal dynamics taking over, b) The transition of the quantum 
walk to the quantum stochastic walk is shown. 



as a quantum stochastic walk, the environment is chang- 
ing the effective graph. Since the quantum stochastic 
walk is more general than the quantum walk, it is there- 
fore universal for quantum computation. Its quantum to 
classical transiton can also be used to examine decoher- 
ence phenomena in quantum computation and classical- 
to-quantum phase transitions. New quantum stochastic 
distributions might suggest new kinds of quantum algo- 
rithms or even new classes of quantum algorithms based 
on this model. 
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